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ABSTRACT 

"O 
O 

We use time series observations from the SOHO and Yohkoh spacecraft to 

\ study solar polar rays. Contrary to our expectations, we find that the rays 

are associated with active regions on the sun and are not features of the polar 

coronal holes. They are extended, hot plasma structures formed in the active 

CN I regions and projected onto the plane of the sky above the polar coronal holes. 

We present new observations and simple projection models that match long-lived 

' polar ray structures seen in limb synoptic maps. Individual projection patterns 

f~ — , ' 

^ ' last for at least 5 solar rotations. 

CO 
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Subject headings: Sun: corona - Sun: UV radiation - Sun: X-rays, gamma rays 
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O ' 1- Introduction 

Polar coronal rays are bright ray-like or arch-like features first seen above the poles of 
the sun during total solar eclipse. They have been described in the literature for at least a 
century (e.g. [Maunder 1901| ) and they are often considered as features of the solar polar 

a' 



coronal holes. Several studies have been conducted to determine the physical properties 
of these high latitude rays ([Fisher and Guhathakurta 19951 , puhathakurta, Fisher and 



Btrong 1996|) . While the classical literature tends to use the terms "ray" and "plume" 



interchangably, the recent studies show that the polar rays have physical properties distinct 
from those of polar plumes. For example, the polar rays are hotter than the background 
polar hole gas (2.6 x 10 6 K vs. 0.7 x 10 6 K; |Guhathakurta, Fisher and Strong 1996| ) while 
the plumes "are cooler than the background by up to (but not more than) 30%" ( Pe Fore"st 



|et al 199"Tp . The reasons for the differences between polar rays and polar plumes have not 
yet been closely examined. Here, we use time series observations made with EIT/SOHO 
and SXT/Yohkoh between 1996 and 1999 to study the morphology and time dependence of 
polar rays. 
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2. Data Analysis 

The observations are based on extreme ultra-violet and soft X-ray full sun disk images 
taken by the EIT and SXT instruments onboard SOHO and Yohkoh, respectively. About 4 
EIT images per day from 1996 through 1998 were used in our study. Each EIT image was 
processed for flat-fielding and filter-correction. The processed images were used to build 
limb synoptic maps. The emissions were extracted around the solar limb 3° wide in polar 
angle and 0.03R Q in altitude. Successive images provide the time domain for this study, 
leading to the limb synoptic maps. Figure 1 shows a limb synoptic map made from more 
than 2000 EIT image pairs 195/171 (FeXII/FeIX,X) taken between January 1996 and June 
1998. Low temperature areas such as the polar holes appear dark while high temperature 
areas such as active regions are bright in Fig.l. To clearly display both southern and 
northern polar holes the range of polar angles (PA) is extended by 90° from 0° — 360° to 
0° — 450° in Figure 1. We used SXT/Yohkoh data in movie form to confirm the distribution 
of the hot plasma associated with active regions moving across the solar disk. 



3. Results 

Figure 1 shows a number of noteworthy features. The coronal polar holes appear as 
dark bands between 150° < PA < 210° and 320° < PA < 390°, extending from day of year 
DOY to about DOY 450 (c.f. |Li et al. 20001 ) . For DOY > 450, the polar holes are filled 



in by a network of multiply crossing rays. Brighter equatorial gas fills the remaining PAs. 
An active region (PA ~ 100° on the east limb and ~ 260° on the west limb) persists over 13 
full rotations from DOY= to DOY= 360. Fingers of emission encroach from low latitudes 
on the polar hole boundary. After DOY= 210, these fingers bridge the polar holes, forming 
a sinusoidal feature that persists for at least 5 rotations. This is dramatically longer than 
the ~ 20 hour lifetime reported for polar plumes ( [Llebaria et al. 1998| ). Sinusoids appear 
simultaneously in the North and South polar holes. At any given time after DOY~450, the 
polar holes are crossed by several overlapping sinusoids making a study of the morphology 
more difficult. In Fig. 2 we enlarge a portion of the limb synoptic map from 210 < DOY 
< 350 days through EIT 284 (FeXV) in which the ray morphology is particularly clear. The 
polar ray links successive appearances of the active regions on the East and West limbs, 
and clearly draws a connection between the polar rays and the active regions. 



To model the polar ray sinusoids, we consider the apparent sky-plane motion of a radial 
coronal structure as it is carried round by the solar rotation. The position angle of such a 
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Fig. 1. — Coronal limb synoptic map at 1.0 R from EIT 195/171 (FeXII/FeIX,X) images in 
the period January 1996 to June 1998. The Day of Year (DOY) is indicated along the lower 
x-axis. The dates are marked along the upper x-axis. The polar angle (PA) is on the y-axis. 
Zero limb angle corresponds to the north pole (NP), while the east limb is at 90° (EL). To 
display the entire southern (180°) and northern (0° and 360°) polar holes, the polar angle is 
extended by 90°. The heliographic latitude is plotted with a solid curve. 
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Fig. 2.— A portion of the synoptic chart from EIT 284 A from DOY= 210 to DOY=350, 
corresponding to July to December 1996. The dashed curve shows Eq. (1) with 9 = 60°, 
the dotted curve corresponds to 6q — 10°. The notations on the right side of the figure have 
the same meaning as those in Fig.l. NP is Northern Pole, EL is East Limb. 
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structure is given by 

tan(6> ) 



ip = tan 1 



cos a 



where #o is the heliocentric latitude and a is the longitude of the feature measured from the 
east limb, 

a = ^(t„ + *) (2) 

Here, P{9q) is the coronal rotation period at latitude 9q, t is time and to is the time when 
the feature crosses the east limb. In Eq. (1) we have neglected the effect of the tilt of 
the solar equator relative to the ecliptic plane. By inspection, we take P{9q) = 27.5 days 
(compatible with independent measurements by [Ingester et al. 19991) . Eq. (1) is plotted on 
Fig. 2 for two sample values 9o = 10° (dotted curve) and 60° (dashed curve). The 9q = 60° 
curve resembles the polar sinusoid rather closely, and suggests that this feature is due to a 
nearly radial coronal ray rooted outside the polar hole but carried across it (in projection) 
by the solar rotation. Models with 9q > 60°, the latitude of the boundary of the polar 
coronal hole at this epoch, do not fit the shapes of the polar sinusoids. The 9q = 10° curve 
fits the polar sinusoid less well but encompasses a broad swath of near-equatorial emission 
that is associated with the active regions. Together, the two curves strongly suggest that 
the polar rays are merely projections of more equatorial gas above the limb of the sun at 
high latitudes. 

A more detailed but still simplistic 3-dimensional corona model was used to further 
test this possibility. Examination of individual images from EIT and SXT suggested the 
following representation. We take the active region to be a source of hot gas, which occupies 
a thick sheet of longitudinal extent A<x The boundaries of the region in the meridional 
plane were assumed to be parabolic with a vertex at the active region. The radial variation 
of the electron density is assumed to follow 



n e (r) = n eO (R )exp 



1.5 x 10 7 ( R 



T e V r 



(3) 



which is consistent with hydrostatic equilibrium. We take T e = 2.5 x 10 6 [K] and 
n e o(R©) = 3 x 10 9 [cm -3 ]. The assumption of hydrostatic equilibrium is not critical to our 
result; any extended distribution of gas that projects above the pole in the plane of the 
sky gives the same qualitative features. The coronal model outside the active region is 
taken from our earlier work ( |Li et al. 2000| ) and consists of separate parts for the equatorial 



regions and polar holes. We integrated along the line of sight (including the effects of the 
tilt of the solar rotation axis) to construct model images of the corona as functions of the 
solar rotation angle. Figure 3 shows a sample model image in which a = 30°. Finally, we 
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constructed limb synoptic maps from the models in the same way as for the data. Examples 
are shown in Figure 4. 

In general, the 3-dimensional models behave exactly as expected from simple 
geometrical considerations (Eqs. 1 and 2). When the heliographic latitude Bq = 0, gas from 
an active region is visible in projection above the polar holes provided it reaches a radial 
distance 

r = , Ra (4) 

y cos 2 #o c °s 2 a + sm #0 

where 9q is latitude and a is longitude as before. A plasma structure on the meridian 
(a = 90°) at 6>o = 10° must extend to r ~ 5.8 R in order to project above the pole, while 
a feature at 9 = 60° needs to rise only to r ~ 1.15 R Q in order to be seen. Therefore, 
the polar rays are most likely the high latitude edges of extended equatorial structures 
seen in projection. We used EIT 171 A (Fe IX,X), 195 A (Fe XII) and 284 A (Fe XV) 
images to estimate the temperatures of the polar rays (e.g. Fig. 1). We find that the 
rays are hotter than the background polar hole gas, consistent with independent results 
( |Guhathakurta, Fisher and Strong 1996|) and with an origin in more nearly equatorial 
regions. The observation that polar rays appear simultaneously in the North and South 
holes (see Fig. 1) is naturally explained in this model if the active region is very extended 
in latitude, consistent with direct images of the active regions when on the solar limb. 

Close inspection of the polar sinusoids in Figures 1 and 2 reveals an additional, more 
subtle long term variation in morphology. The rising (PA increasing with time) and falling 
(PA decreasing with time) parts of the polar sinusoids have unequal brightness. In the 
interval 210 < DOY < 350, the rising branch is systematically brighter than the falling 
branch, for both northern and southern features. The sense of this asymmetry reverses 
in the period 350 < DOY < 530 (Fig. 1). In fact, Fig.l shows that the asymmetry 
alternates in phase with the heliographic latitude (B Q ), suggesting that projection effects 
are responsible. 

Our 3D coronal model, which includes Bq, shows exactly this behavior (see Figure 
4). When the north pole is tilted towards the Earth (Figure 4a), the model polar rays are 
brighter on the rising branch than on the falling branch. The opposite pattern is reproduced 
in Figure 4b when the north pole is tilted away from the Earth, providing a realistic match 
to the data. When the north pole is tilted towards the Earth, the hot active region gas 
must rise higher on the far side of the sun in the southern hemisphere (and on the near 
side in the north) to meet the line of sight, and so appears fainter than front-side gas on 
the line of sight at lower altitudes. Therefore, the bright rising segments near DOY ~ 242 
(see Fig. 2, from August 23 to September 4, 1996) identify an active region on the visible 



Fig. 3. — Sample model image in which the active region is located at 30° longitude and the 
north pole of the sun is tilted towards Earth by 7.°25 (i.e. B = +7.°25). Image has north 
to the top, East to the left, and the grid spacing is 10°. 
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Fig. 4. — Simulated limb synoptic map made from the model of Figure 3. A rotational 
resolution of 1° was used. About 5 solar rotations are included here. The panels show a) 
B = +7.°25 and b) B = — 7.°25. Note the reversal in the sense of the asymmetry of the 
brightness of polar rays. 
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hemisphere of the sun, as is seen. The periodically reversing brightness asymmetry is a 
product of projection effects in equatorial plasma. 

As the solar cycle approaches its activity maximum, the low temperature gas in the 
polar holes is progressively obscured from view by hot equatorial gas projected to high 
latitudes. This is well seen in Fig.l where, after DOY ~ 550, both polar holes are covered by 
a worm-like network of bright polar rays. The basic morphology of individual rays appears 
unchanged in the 2| years of this study: only the number of polar rays increases towards 
solar maximum. Some previous reports have claimed that the polar holes vary in size and 
strength through the solar cycle. Our present results suggest instead that the morphology 
and visibility of the polar holes are largely controlled by the number and strength of active 
regions far outside the holes and by the associated high altitude gas that projects across 
the holes as polar rays. In future work, we hope to explore this possibility. 



4. Summary 

1) Limb synoptic maps of the sun in the period 1996-1998 show long-lived, periodic 
features in the inner solar corona. 

2) Polar rays appear clearly in such maps with a morphology that shows them to 
be hot gas structures from active regions projected above the polar holes. Unlike polar 
plumes, which are clearly rooted in the polar coronal holes, the polar rays are not physically 
associated with the holes. 

3) The correlated appearance of polar rays above both north and south holes shows 
that these features are caused by equatorial plasma structures with a latitudinal extent 
comparable to the solar diameter. 

4) Annually reversing asymmetries in the polar ray brightness on limb synoptic maps 
are caused by projection effects due to the varying heliographic latitude. 

5) Individual polar rays are long lived (more than 5 rotations in this study) and 
increase in number as the solar activity cycle approaches maximum. The longevity of the 
polar rays is a reflection of the longevity of active areas that are large enough to be seen in 
projection above the solar poles. 

6) Variations in the visibility of the coronal polar holes are largely controlled by 
projection of hot gas from more equatorial regions. Polar plumes are features best seen at 
solar minimum. Polar rays are most prominent near solar maximum. 
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